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We explore several strategies of eliminating (or nulling) the small-scale information in weak lens- 
ing convergence power spectrum measurements in order to protect against undesirable effects, for 
example the effects of baryonic cooling and pressure forces on the distribution of large-scale struc- 
tures. We selectively throw out the small-scale information in the convergence power spectrum that 
is most sensitive to the unwanted bias, while trying to retain most of the sensitivity to cosmological 
parameters. The strategies are effective in the difficult but realistic situations when we are able 
to guess the form of the contaminating effect only approximately. However, we also find that the 
simplest scheme of simply not using information from the largest multipoles works about as well as 
the proposed techniques in most, although not all, realistic cases. We advocate further exploration 
of nulling techniques and believe that they will find important applications in the weak lensing data 
mining. 



I. INTRODUCTION 

Progress in cosmology in the last decade has been dra- 
matic, and there is reason to believe that we will make 
further advances on key questions in the coming years. 
Key to this progress is the ability to confront precise 
data with highly accurate theoretical predictions. Grav- 
itational lensing illustrates both the potentials and diffi- 
culties of this era of cosmological research. In principle, 
lensing can constrain the properties of the dark energy 
causing the accelerated expansion of the universe and 
constrain both the cosmic geometry and thegrowth of 
large scale structure (e.g. 0, II 13, 0, IE S 01a H ) ■ But 
one of the main challenges for future weak lensing sur- 
veys will be controlling the theoretical systematic errors 
involved in predicting the lensing signal at small, nonlin- 
ear scales of <C 1 degree. These errors include numeri- 
cal artifacts in c omp uting the non-linear power spectrum 
(see e.g. 0,0, and complex baryonic processes 
that are difficult to model accurately [La, fl5l | and make 
semi-analytic predictions of nonlinear power unreliable. 

Many of the difficulties in making theoretical predic- 
tions come in at small physical scales. In this paper we 
explore one approach to protect against theoretical bi- 
ases at small scales. We suggest selectively throwing out 
the small-scale information and "nulling" the bias while 
at the same time retaining useful cosmological informa- 
tion. What makes us optimistic in this regard is that the 
information contained in the weak lensing convergence 
estimated in different redshift slices is strongly overlap- 
ping, essentially because of the width of the lensing ker- 
nel. Therefore, dropping a reasonably small subset of the 
tomographic information leads to small degradations in 
overall cosmological constraints. This has been used by 
Takada & White 0] to propose dropping the autospec- 
tra in order to protect against unwanted biases due to in- 
trinsic alignments of galaxies. We take this general idea 
further and propose identifying specific linear combina- 
tion of the cross-power spectra that are most sensitive 
to a given nonlinear effect, and then dropping them in 



the analysis in order to null out the effect while ideally 
preserving most of the sensitivity to cosmological param- 
eters. 

In the next Section we describe the fiducial assump- 
tions about the survey and small-scale bias. In Sec. IIIII 
we describe and explore several nulling tomography tech- 
niques. We discuss the results and comment on future 
prospects in Sec. IVII 

II. METHODOLOGY AND SMALL-SCALE BIAS 
MODEL 

Assume a weak lensing survey with the ability to di- 
vide source galaxies in N s redshift bins. Let PMi) be 
the convergence cross-power spectrum in ith and jth to- 
mographic bin at a fixed multipole I — for definitions 
and details, see e.g. Ref. 0- The observed convergence 
power is 

C? j {l)=P%{l) + 8 ij { ^ (i) 

where (7j„ t ) 1 ^ 2 is the rms intrinsic shear in each com- 
ponent which we assume to be equal to 0.22, and fii is 
the average number of galaxies in the ith redshift bin per 
steradian. 

For definiteness we assume a SNAP-type survey 0] 
covering 1000 sq. dcg. with the galaxy distribution of 
the form n(z) cx z 2 exp(— z/zq) that peaks at 2zq = 1.0. 
We assume 100 usable galaxies per arcmin 2 . We use in- 
formation from the wide range of scales corresponding 
to multipoles 100 < I < 10, 000. We have a set of six 
cosmological parameters: energy density and equation of 
state of dark energy Ode and w, spectral index n, matter 
and baryon physical densities f^M^ 2 and fish 2 , and the 
amplitude of mass fluctuations <7g. Throughout we as- 
sume a flat universe. For our fiducial results, we let these 
six parameters vary without any priors since the fiducial 
survey alone is powerful enough to determine the cosmo- 
logical parameters to a good accuracy. The cosmological 
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constraints can be computed using the standard Fisher 
matrix formalism. 

To be definitive, let us now consider systematic bias in 
the matter power spectrum P(k) of the form 



P(k )Z )=P(k,z) [l+o/(fc)] 



(2) 



where P(fc, z) is the true, unbiased power spectrum, and 
a is the coefficient that we use to "turn on" the bias 
f(k) whose functional form is, we assume for a moment, 
precisely known. [Note that realistic biases f(k) may 
be redshift dependent as well, perhaps increasing at late 
times during the nonlinear structure formation, but we 
ignore this issue since it is unimportant in the follow- 
ing discussion.] In order to use weak lensing to measure 
the cosmological parameters, we are forced to fit for the 
parameter a as well, otherwise our results will be biased. 
Our goal here is to minimize the sensitivity of the data on 
parameter a, since the results would thus become less sen- 
sitive to the particular bias in the power spectrum from 
Eq. (J2J). At the same time, we hope not to significantly 
weaken the sensitivity to other cosmological parameters, 
such as the equation of state of dark energy w. 

For the bias in the matter power spectrum, our fiducial 
model is as in Eq. J2J) with a = 1.0 and f(k) = (fc/fc*)". 
This model phenomenologically describes the baryonic 
effects if we further set a ~ 3 and fc* ~ 10/iMpc -1 0, 
IbH ] and these are the values we use, although we explored 
a range of other values and found similar results. Finally, 
to compute the bias in the cosmological parameters due 
to the small-scale bias, we follow the standard formalism 
that uses the Fisher matrix described in e.g. Ref. [l3| . 



covariance matrix of those combinations is 

(uiUj) = QikQjlCkl = q (i)T C q« (5) 
k,l 

where the column vector qW is defined as the i th row 
of the matrix Q, and the matrix C has elements = 
{l(i,m),kl(t,m),i) ■ Say we only had a measurement of a 
single Ui = u. The likelihood for u (in the Fisher matrix 
approximation) is gaussian 



C 



1 



: exp 



2R 



(6) 



where (u 2 ) = R = q T C q. and we suppressed the index 
i on q as well. Then we want to determine the q by 
maximizing the error in a, or minimizing (after a bit of 
algebra) 



OR 
da 



(7) 



We now need to solve the following problem: find the col- 
umn vector q that maximizes F aa in Eq. JJJ, given the 
N s x N s covariance matrix C. Fortunately, this problem 
is well known and can be solved using standard tech- 
niques (e.g. 0,0|). The solution is to consider 



da 



(L 



-1\T 



(L T c 



A(L T q) 



(8) 



where L is a lower triangular matrix satisfying C = LL T . 
Eq. JSJ is an eigenvalue-eigenvector relation for L T q 
which we can solve to obtain q. It can easily be shown 
that the error in a from measurement of any single com- 
bination v., is v2/|Ai|, and that (iiiUj) — Sjj. Therefore, 
we get the whole spectrum of uncorrelated combinations 
of shear ordered by their sensitivity to the small-scale 
parameter a. Again, recall that we need to repeat the 
same procedure at every multipole £ (or, multipole bin 
centered at £). 

We find that the Fisher information on the parameter 
a is heavily dominated by one (or a few) particular com- 
binations of shears. For example, for the multipole bin 
centered at £ = 1100 and 4-bin tomography, the eigenval- 
ues A are: 7.0 • 10~ 6 , 4.2 • 10~ 5 , 3.3 • 10~ 4 , and 1.3 • 10~ 2 . 
Clearly, by far the most information about the nuisance 
parameter a is carried in the eigenvector corresponding 
to the fourth eigenvalue. This eigenvector is (after ad- 
justing the unimportant overall normalization for clarity) 
7x — 5.372 + IO.273 + 5.874. Therefore, this combination 
of shears at £ = 1100 is most sensitive to the parameter 
a and would be the one to drop. 

Figure ^] shows dC^/dhxk and d(uiUj)/d\nk at mul- 
tipole £ = 1100, now for a 10-bin tomography, where 
are the original cross-power spectra and (uiUj) are the 
covariances of the linear combinations of shears Ui de- 
fined in Eq. It is clear that the covariances of the ui 
have a more differentiated wavenumber dependence than 
where hereafter we drop the subscript (£, to) from i^. The the original C" L However, because of the width of the 



III. 1-POINT NULLING 

Let us first consider using N s specific linear combi- 
nations of shears in individual redshift bins, specifically 
Ylf=i Qij'Yjifyi where the coefficients Qij are indepen- 
dent of the angular position of the galaxy. In other words, 
each linear combination gives equal weight to all galaxies 
in the jth redshift bin. Our goal is to find specific linear 
combinations that are most sensitive to the nuisance pa- 
rameter a and simply throw them out. The nice feature 
of this problem is that it is tractable analytically. 

Let us fix the multipole £ and define 7(i jrra ),i = C(^,m)ii 
where ar£ t m) t i is the spherical harmonic coefficient of the 
shear map in redshift bin i. The only things we really 
need to know about 7^^)^ are its mean and variance 



Now consider the linear combinations 
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FIG. 1: Derivatives dC^/dlnfc (left panel) and d(uiUj)/d\nk (right panel) at multipole £ = 1100. Here C£ are the original 
cross-power spectra and (itiUj) are the covariance of the new combinations of shear, denned in Eq. ©. Note the coverage in k 
for the new eigenmodes; the ones that have most of their weight at high wavenumbers are most sensitive to the parameter a. 
The overall normalization in each panel is arbitrary, but the relative size of the modes is not. 



gravitational lensing kernel, the weights d(uiUj) fdliik, 
shown in right panel, are still strongly overlapping. 

First, consider the fiducial case where we keep all infor- 
mation in the survey, but add a single new parameter a 
describing the small-scale effect. If the form of the small- 
scale effect is known exactly, we find that the additional 
degeneracy introduced by the new parameter a is negligi- 
ble, and marginalized errors in cosmological parameters 
increase by only a few percent relative to the case when a 
is fixed. This is not surprising, as the parameter a enters 
the observables in a different way than the cosmological 
parameters, and the fiducial survey is powerful enough to 
determine them all without substantial loss of accuracy. 

Let us now drop the combination Ui that is the most 
sensitive to parameter a. We find that the marginalized 
error on the parameter a increases by several orders of 
magnitude while increasing the cosmological parameter 
accuracies only by several tens of percent. While this 
sounds like fantastic news, we should remember that we 
optimistically assumed that we were able to exactly guess 
— and parametrize — the form of the small-scale effect 
f(k). In reality, we will be able to guess f(k) only approx- 
imately, and the difference between the true and guessed 
f(k) will lead to biases in the measured cosmological pa- 
rameters. The real figure of merit is the ratio of the bias 
due to the incorrectly guessed f(k) and the fiducial ac- 
curacy in any given cosmological parameter. Clearly, our 
goal should be to strike balance between the systematic 
bias and the statistical precision, minimizing the former 
while not increasing the latter by more than a few tens 
of percent. 

Table [I] shows the errors and biases in w when be- 
tween zero and two combinations of shear, out of ten 
total, were dropped, and alternatively, when multipoles 
(■ > 4iax were dropped. While our 1-point nulling works 
impressively well, the cut in multipole space is as, and 



perhaps even more, impressive. For example, to achieve 
the same statistical error as in 1-point nulling with one 
combination of shears dropped (error > 100% bigger than 
the fiducial one) we can alternatively cut all multipoles at 
i > 1000, but then the bias with cutting in I is more than 
two times smaller than the bias with nulling! We have 
checked that these results are qualitatively unchanged 
with a different choice of both the actual and the guessed 
bias. 



IV. 2-POINT NULLING 

Inspection of Table |U shows that one problem with 1- 
point nulling is that it does not have enough resolution in 
choosing how much of the small-scale structure to null. 
For example, dropping the single combination of shears 
that is most sensitive to the parameter a increases the 
error in w by about 150%, which is too much to tolerate 
regardless of the success in nulling out the unwanted bias 
in P(k,z). The lack of ability to perform nulling more 
gradually is not too surprising, as there are only N s linear 
combinations to choose from. 

Ideally we would like to have a better resolution in 
the eigenvalues Xi and thus obtain a better leverage in 
controlling the cosmological parameter degradations as 
a function of the nulling efficiency. One way to remedy 
that is to form linear combinations of the convergence 
power spectra C*-, as there are N 
them 



N S (N S + l)/2 of 



N 



(9) 



where Qij are coefficients (different from Qij in the pre- 
vious section) and we denote a pair of redshift bins by 
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1-point Nulling 


Cutting above Imax 


Skipped 


a(w) 


5{w) 


f 

*-max 


cr(w) 


5{w) 





0.054 


0.085 


10000 


0.054 


2.630 


1 


0.130 


0.023 


3000 


0.073 


0.101 


2 


0.291 


0.003 


1000 


0.117 


0.010 



TABLE I: Nulling out shear combinations that are most sen- 
sitive to a vs. simply cutting measurements at some i?max for a 
10-bin tomography case. The nulled SP(k) was assumed to go 
as k 3 (with k* = lOh Mpc -1 ), the actual one as k 35 . While 
nulling in principle works well, it leads to degradation in sta- 
tistical errors that are about 150% or larger, and therefore 
too large to tolerate. On the other hand, simply throwing out 
power beyond some £ max is surprisingly efficient; for example, 
restricting to ^ ma x = 1000 leads also to a similar > 100% in- 
crease in ct(w), but with the bias in w that is smaller than 
with 1-point nulling. More importantly, cutting in £ can be 
done gradually, while the 1-point nulling necessarily leads to 
jumpy changes in the statistical error and bias. 



a single subscript; here i and j take values from 1 to 
N. The goal is to find the linear combination(s) Vi that 
produce the largest Fisher matrix element F aa . While 
the 1-point procedure can be repeated, the optimization 
problem, sadly, cannot be solved analytically as before 
because (Cj) ^ 0, and we needed that assumption for 
the expression for F aa , Eq. Q, to take its simple form. 
Therefore, we choose to find the optimal combinations 
of the convergence power spectra by brute force. Us- 
ing Powell's minimization method |20|. we find the lin- 
ear combination v\ that has maximal F aa . Using the 
Gram-Schmidt algorithm, we then project to the (N — 1) 
dimensional subspace, perform the maximization again, 
and find the second most sensitive combination, V2- And 
so on, until we find all combinations of the power spec- 
tra ordered by their sensitivity to the small-scale effect. 
As with the 1-point function, we apply this algorithm 
at every £ separately, where £ is the center of the corre- 
sponding band power window in multipole space. 



We indeed find that the resolution in nulling is greater 
with 2-pt nulling than with 1-point; for example, throw- 
ing out the most sensitive combination increases a(w) by 
about 40%, and not 150%. However, we found two prob- 
lems with 2-point nulling. First, as before we find that 
simple cutting in £ works about as well or better. In ad- 
dition, we find that the biases in cosmological parameters 
are substantially suppressed only if the form of the bias 
is nearly precisely known, and not in more realistic situa- 
tions when the bias is guessed incorrectly. These conclu- 
sions are unchanged even after trying different forms of 
f(k), looking at different cosmological parameters, differ- 
ent values of fc* , and different numbers of survey redshift 
bins. Therefore, unless our 2-point procedure is substan- 
tially improved in some way, 1-point nulling is the more 
effective approach. 



V. CUTTING IN WAVENUMBER k 

Another simple strategy is to simply throw out contri- 
butions to the band powers coming from wavenumbers 
greater than some (comoving) cutoff fc cut . Assume that 
the Limber's approximation integral over the large-scale 
structures is discretized, with N p (rather than infinitely 
many) lens planes. Then, assigning a single subscript i 
for any of the N cross-power spectra, we have 



3=1 



(10) 



where W are suitably defined weights that depend on 
the lensing geometry and the distribution of galaxies and 
P(kj, Zj) is the matter power spectrum evaluated at some 
wavenumber kj — £/r(zj) and redshift Zj of the lens plane 
in question. We now form a, N x N p matrix with the 
entry equal to WijP(kj, Zj); we have N = 10 (cor- 
responding to N s = 4 redshift bins) and N p = 59. Then 
we simply transform this matrix to the lower triangu- 
lar form. The rows now represent linear combinations 
of the convergence power spectra, and new matrix en- 
tries are their weights. Of course, all N rows still contain 
the same information as before, as we simply performed 
a linear operation on the C K s. However, by consider- 
ing the first M rows only (M < N), we can avoid using 
any information coming from wavenumbers greater than 
the wavenumber on (N p — M)th lens plane. As before, 
we repeat the procedure at each multipole band centered 
at I. For simplicity we hold fc cu t independent of £, and 
then higher I implies higher kj for a fixed lens plane — 
therefore, more combinations will be cut at higher t, un- 
til some multipole i ~ fccutf (^max) where no information 
at fc > fc cut will be left. Note too that, unlike the 1 and 
2-point nulling, cutting in fc does not require a guess of 
the offending small-scale effect as all information beyond 
fccuti "good" and "bad", is thrown out. 

To explore the efficacy of this type of nulling, we per- 
form it for a variety of values fc cu t, then compare to 
the "gold standard" of dropping multipole power spec- 
tra above a fixed € max . Overall, we find that cutting in 
k works better than cutting in £, as expected since our 
small-scale contamination is defined in k. However, the 
two methods are essentially equally powerful unless the 
small-scale effect enters at very large scales (i.e. if fc* is 
small), and in that case the fc-cut can be significantly 
more powerful; see Fig. [21 This is easy to understand as, 
from the relation fc ~ £/r(z), we see that even if stringent 
cutting in £ is performed, some amount of high fc modes, 
corresponding to low redshift structures, will be left in 
the data — and these high fc modes can be extremely 
damaging if fc* is small since the bias in the observables 
goes as a power law in fc/fc*. It is also interesting to 
note that cutting in fc is equivalent to removing the low- 
redshift structures from the data — in other words, in- 
stead of cutting above some fc cu t , one could equivalently 
choose not to use lens galaxies with (photometric) red- 
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FIG. 2: Comparison of the cutting in wavenumber k vs. the cutting in multipole I. Left panel: evolution of error vs. the 
(absolute) bias in J7de with increased aggressiveness in cutting in fc-space (solid black curve) or ^-space (dashed blue curve). 
The cut is increasingly more severe going from bottom right to top left. Here we use a 4-bin tomography case with the small- 
scale effects going as (fc/fc*) 3 with with fc* = 10/iMpc -1 . For orientation, we show one value each of ^ ma x and fc cu t on the curves 
along which they vary; these are roughly the points where the bias is reduced to a fraction (~ 0.5) of the statistical error, and 
in this case the total error with the fc-cut is about 50% smaller than that with the ^-cut. Right panel: statistical error in Qde 
after cutting as a function of the scale at which the small-scale bias enters, fc» , for the k and i cutting. With each method and 
for each value of fc*, just enough cutting is performed so that the remaining bias in Que is equal to one half of its statistical 
error. Note that the two methods are comparable if the bias enters at small scales (high fe*), but if the bias affects much larger 
scales the fc-cutting is clearly more effective as it leads to a smaller final error in J7de- 



shifts less than z cu t where r(z cu t) — £/k cu t- Which of 
these two possibilities is more feasible in practice is an 
interesting question, but outside the scope of this work. 

Current estimates indicate that the effective scale of 
various small-scale effects is anywhere from a few to a few 
tens of hMpc- 1 (e.g. [H Q If fe* is closer to the 
lower end of this range for any of these effects, cutting 
out the high wavenumbers (or low-redshift structures) 
may prove extremely beneficial. 



VI. DISCUSSION AND CONCLUSIONS 

We have described and explored three different strate- 
gies for nulling out small-scale information in weak lens- 
ing power spectrum measurements, with the idea of pro- 
tecting against uncertain, approximately known contam- 
ination while preserving most of the sensitivity to cosmo- 
logical parameters. 

The 1-point nulling tomography uses select linear com- 
binations of shear in different redshift bins so as to max- 
imally reduce the survey's sensitivity to small-scale bias. 
The algorithm for achieving this is well known and largely 
analytic, and it produces the desired linear combinations 
ordered by how sensitive they are to the parameter a 
describing the small-scale bias. Figure ^ shows that, 
as expected, the combination that is most sensitive to 
the small-scale bias has weight at smaller scales than the 
other combinations. Since in a realistic situation we will 
only have approximate understanding of the true bias, 



we test the performance of 1-point nulling by assuming 
the fiducial bias is a power law in wavenumber, while 
the fit bias, normalized by the parameter a, is a different 
power law. We find that even as just one linear combina- 
tion of shear is dropped from the analysis (the one most 
sensitive to the parameter a) the biases in cosmological 
parameters are sharply reduced. However, in that case 
the increase in cosmological parameter errors is about 
150%. Whether that is too large to tolerate depends on 
the fiducial power of the survey and on the deleterious 
effect of other, unrelated systematic errors. 

The problem of large initial degradation can in prin- 
ciple be ameliorated by 2-point nulling, where we form 
linear combinations of the cross power spectra (rather 
than shear). Because the number of cross power spectra 
is of order the number of redshift bins squared, we have 
many more observables at our disposal and can achieve 
a finer resolution in the amount of nulling performed, 
and hence in cosmological parameter error degradation. 
Unfortunately 2-point nulling does not have a nice ana- 
lytic solution and linear combinations of power spectra 
that are most sensitive to the small-scale bias need to be 
found by brute force. We find that the method works 
well in cases when the functional form of the bias can 
be guessed accurately. If the exact form for the bias is 
guessed incorrectly, however, 2-point nulling produces bi- 
ases in cosmological parameters that do not decrease as 
the most sensitive combinations are dropped. We con- 
clude that 1-point nulling is more effective than the 2- 
point approach, and this is unlikely to change unless 
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a substantial improvement to our 2-point procedure is 
found in the future. 

While both 1 and 2-point nulling are effective, we found 
that they are not significantly better than the simplest 
strategy of dropping the highest multipoles of the con- 
vergence power spectra. [This may not be too surpris- 
ing in retrospect: decreasing £ max from 10000 to 3000, 
for example, increases the cosmological parameter error 
bars by only ~ 30% while having a tremendous impact 
in removing a variety of small-scale contaminations - 
in other words, the price to pay with dropping the high- 
est multipoles is relatively small, and this procedure is 
already reasonably effective.] This motivated us to ex- 
plore yet another, but very different, method for remov- 
ing the small-scale biases. We proposed a method of 
removing all information above some wavenumber A; cut 
or, roughly equivalently, removing structures below some 
redshift from the data. This approach is different from 
the previous two in that the functional form of the bias 
need not be known, and all information above fc cut , useful 
or not, is thrown out. This works the best of all methods, 
leading to cosmological parameter errors that are 10-50% 
percent smaller than those with the I space cut, and by 
up to a factor of two if the small-scale effect enters at suf- 
ficiently large scales. Given that weak lensing shear data 
is likely to be biased by large nearby structures (e.g. a 
large cluster of galaxies at low redshift), some form of 
cutting in wavenumber will surely be beneficial to pro- 
tect against biases. 

All of the aforementioned techniques require rough 
knowledge of the cosmological model so that the fiducial 
convergence power spectra (i.e. those not including the 
small-scale bias) are known to a good accuracy, and can 
be manipulated to null out the unwanted effects. This 
should not be a problem in the future, as the cosmolog- 
ical parameters arc already reasonably well determined 
with the combination of cosmic microwave background, 
type la supernovae and large-scale structure constraints. 



However, systematic errors in weak lensing measurements 
are a bigger concern. While the study in this paper was 
precisely concerned with removing biases due to a class of 
systematics - uncertain theoretical predictions on small 
scales - there will be other systematics in weak lensing 
measurements that can roughly be divided into redshift 
errors and additive and multiplicative errors in measure- 
ments of shear | 21| . Those errors are likely to depend on 
numerous factors (e.g. observational season, galaxy mor- 
phology, etc.) and are uncertain at this time. We do 
not believe that the systematics will significantly affect 
the nulling techniques described in this paper since the 
systematics, provided they are small, should be a lower- 
order contribution to our observables (which are linear or 
quadratic in shear). Nevertheless, just as with all other 
applications of weak lensing in the next generation ex- 
periments, incorporating the systematics into the nulling 
tomography will be an important and challenging task. 

In conclusion, we strongly believe that the strategies 
that we described will find application in weak lensing 
data mining. We only considered the convergence power 
spectrum information here. Other aspects of weak lens- 
ing, such as mapping of dark matter, depend more criti- 
cally on the small-scale information, and it might be that 
nulling-type techniques will reach their full potential in 
precisely those applications. It is likely that these inter- 
esting possibilities will be considered in the future. 
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